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Abstract 

The shape invariance condition is the integrability condition in supersymmetric quantum mechanics (SUSYQM) . It is a 
i—{ • difference-differential equation connecting the superpotential W and its derivative at two different values of parameters. We 
■ show that this difference equation is equivalent to a non-linear partial differential equation whose solutions are translational 
00 ■ shape invariant superpotentials. In lieu of trial and error, this method provides the first ab initio technique for generating 
shape invariant superpotentials. 
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Supersymmetric quantum mechanics (SUSYQM) [T] extends Dirac's factorization methoqfl for the harmonic os- 

Q_i; 

cillator to a large number of other potentials, whose energy eigenvalues can then be obtained algebraically. The 
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SUSYQM extension consists of introducing a superpotential W(x, a) that generates two partner Hamiltonians, 
both with the same energy eigenvaluejf). These partner Hamiltonians are given by 

H T = A ± A* = (t h ^ + W(x, a)j(± h ^ + W(x, a] 

*2 d 2 TTr7 , , . dW(x 7 a) _2 d 2 Tr . . . . 

= -Jt 2 -^+W 2 (x,a)Th d -^ = - fl V + T(l,a) ' () 

where partner potentials V±(x, a) are related to superpotential W by W 2 (x, a)±h dW ^°^ ■ (We have set 2m = 1.) 
One of the best known examples is the superpotential W(x, a) — —a cot a; with a > 0, and < x < tt. The partner 
potentials generated by this W(x, a) are: V±(x, a) = W 2 (x, a)±h dW ^' a ^ = a(a±Ti ) cosec 2 x — a 2 . For the special 
case of a = h, the potential V-(x,h) is a constant function: — a 2 = — Ti 2 ; i.e., just an infinite one-dimensional 
square well potential whose bottom is set to — h 2 , while its partner potential V+{x, K) is given by h 2 (2cosec 2 x— 1). 
Thus, the much less familiar cosec 2 potential has the same spectrum as the well-known infinite well (except for 
the infinite well's 'extra' ground state). 

However, for any pair of partner potentials, one would still have to know the spectrum of one of the partners 
to obtain the other. The critical advance in SUSYQM was the discovery [3] and subsequent rediscovery of shape 
invariance (SI) 4 . A system is called shape invariant if partner potentials could be shown to differ only by the 
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3 Dirac however, credits Fock with that discovery: P. A. M. Dirac, Communications of the Dublin Institute of Advanced Study, Al, pp. 5-7 
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value of a parameter and an additive constant: 

V+(x,a ) = V-{x,ai) + R(a ) , (2) 

in which case each spectrum can be generated without reference to its partner. 
We can write, more symmetrically 

V + (x,a ) + g(a ) = V-(x, a-s) + #(ai) , 

where i?(ao) = g{a\) — g(ao)- That is, 

w 2 ( Xl a ) + h + 9 (ao) = w 2 {x, ai )-n + fl ( Bl ) (3) 

The shape invariance condition in the form of this difference-differential equation has also been studied for 
dynamical systems, where it is known as the infinite-dimensional dressing [6j. In the case of our example above, 
W(x, a) = — a cot a;, the partner potentials are related by 

V+(x,a)+ a 2 = V-(x,a + h) + (a + h) 2 . (4) 

Thus they are shape invariant: a shift of parameter ao = a to a± = a + h <?(ao) = a 2 ,, and an additional constant 
(a + h) 2 — a 2 is all that distinguishes the potentials V+(x, a) and V-(x, a + h) from each other. But we could just 
as well have chosen ao = a + fi and a\ = a + 2fi, and written 

V+(x,a + h) + (a + h) 2 = V-{x,a + 2h) + (a + 2hf . (5) 

By induction, we can generate a family of partner potentials related by a n = a n -\ + h. This form of shape 
invariance is called "additive" or "translational" [1] . This will be the case for all the shape invariant superpotentials 
we shall consider in this paper. (There are other forms of shape invariance, but we shall not consider them here.) 
We emphasize that the shape-invariance condition (Eq. ([3])) is non-linear in W; therefore one cannot simply add 
arbitrary constants to W and maintain shape-invariance. This extremely strong constraint restricts the known 
shape invariant superpotentials to a very small number. 

For shape invariant superpotentials, SUSYQM yields [TJ energy eigenvalues 

Et\a ) = J^R(<H) = g(a n ) - g(a ). 

Thus, for the infinite well spectrum, with its bottom at —a 2 = — fi 2 , 

Et\a ) = g(a n ) - g(a ) = (n + I) 2 h 2 - h 2 3 

Note that Eq. ([3]) a difference-differential equation; that is, it relates the square of the superpotential W and its 
spatial derivative computed at two different parameter values: (x,ao = a) and (x, a\ = a + fi). Our objective 
is to show that this integrability condition can be expressed as a non-linear partial differential equation that is 



5 N.B. Unlike the standard convention in quantum mechanics texts, n begins at 0, the well width is 7r, and the well bottom is at —% . This 
results into a ground state energy ' = 0, a requirement of the SUSYQM formalism. 
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local in nature; i.e., all terms in the equation can be computed at the same point (x, a). This has the obvious 
advantage of mathematical familiarity (at least to physicists). It provides a systematic method for finding the 
superpotential, and thus g(a n ) and from it the energy spectrum of translational shape invariant superpotentials. 
We shall show that all known translational shape invariant superpotentials [1] are its solutions. 

As in the above example we will only consider translational shape invariance, where the parameters are related 
by fii = ao + h- For some superpotentials, such as Morse, successive parameters see a decrease by fi instead of 
an increase. Such superpotentials are also included in these discussions. For the Morse superpotential W(x,A) = 
A — Be~ x , if we insist that A — > A — Ti at each step of change of parameters, A of Morse will be related to the 
a parameter by A — —a. Note that we have made no assumptions about Ti: for now, it is simply a part of the 
parameter defining W. Since this equation is valid for any value of h, we will treat x and a — Ti on equal footing 
as independent variables. Equation ([3]) can now be written as: 

W*( X ,a-h) + h dWix ' a - h) +g(a-h) = W*(x,a) - fc^<M + g{a) ( 6 ) 
ox ox 

where <ij = a and = a — %. If we now differentiate Eq. © with respect to h, noting that W 2 (x, a) and g(a) 
are independent of H, we obtain 

dWj^a-h) dW(x,a-h) 3 dW(x,a-h) dg{a~h) dW(x,a) 
2W(x,a-h) m + ^ +h- + = %r- (7) 



Now because of translational shape invariance, the dependence of W on a and h is through the linear combination 

(x,a — h) dW(x,a 

dh da 



a — h; therefore, the derivatives of W with respect to a and h are related by: dw ( x ' a h ^ — — dw ( x ' a n > 



Then we may rewrite Eq. ([7]): 

2W{x a fy dW ( x > a ~ % ) + dw ( x > a - h ) n ddW{x,a~K) dg{a - h) _ dW{x,a) 
da dx da dx da dx 

But this must be true for any value of fi. If we then set Ti = 0, shape invariance implies the following condition 
for the superpotential W(x, a): 



o w/ ^ f dW(x,a) \ dW(x,a) dg(a) 



(8) 



This is the desired partial differential equation which is equivalent to the original difference-differential equation 
for shape invariance. Thus, we have shown that all shape invariant superpotentials with a^+i = ai + h are solutions 
of the above non-linear partial differential equation. Before we develop a method to determine solutions of this 
complex equation; i.e., to find shape invariant superpotentials, we will as an example show that two known shape 
invariant superpotentials are indeed solutions of this differential equation, and we shall obtain their spectra. It 
can be easily verified that all members of the list given in Ref. [TJ [5] are solution of this equation. 
Morse potential: The superpotential is given by 

W(x,A) = A-Be~ x (9) 
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Substituting this superpotential in Eq. l[8"]). we get: 




whence, A = — ^ q A ' . As stated earlier, we have eto = —A and a n = ao + n% = —A + nh, with g(A) = —A 2 . 
This leads to E n = A 2 - (A- nh) 2 . 



Scarf I potential: The superpotential is given by 



W(x,A) = ,4 tan x + Bsec x A>B, 



(11) 



where A is the parameter that changes. Substituting this superpotential into Eq. ([5]), we get: 



2 (A tan x + Bsec x) (tan x) 



2 (A sec 2 x + B sec x tan + 



dg(A) 







<9A 



3<7(A) 



- 2 A = 



<9A 



hence, <?(j4) = A 2 . Thus in our notation do = A and a„ = do + nh = A + nh, and E n = g(a n ) — g(aa) — 



One can check that all examples listed in Ref. [II [2] satisfy the partial differential equation given in Eq. ([8]). For 
a given shape-invariant superpotential, Eq. ([S]) determines ,g(a), and hence its eigenspectrum can be obtained 
using E n = g(a n ) - g{a Q ). 

However, our true goal is the opposite. Until now, finding shape invariant superpotentials was left to excellent 
intuitions or trial and error. We develop for the first time an ab initio systematic method for these superpotentials. 
The remainder of this paper shows how this can be done. 

Since all translation shape invariant superpotentials are necessarily solutions of Eq. ([5]), finding solutions of 
this equation will help us in our quest for shape invariant superpotentials. For this, we have developed a method, 
similar to the separation of variables, to solve this non-linear partial differential equation. Using our method, we 
are able to derive all the known translational shape invariant superpotentials listed in Ref. [TJ [5] . The method 
also provides a check as to whether there are additional, as-yet-undiscovered superpotentials. 

At this point, it is tempting to solve Eq. ([8]) by assuming various forms for the term ■ However, instead 
of that rather random process, we choose to try a factorizable ansatz 

JV 

W(x,a)=Y,Ma)X l (x) (12) 

i=l 

where Ai(a) and Xi(x) are assumed to be respectively functions of a and x. Thus, to determine W(x, a), we will 
need to find 2A^ functions Ai(a) and Xi(x). Eq. ([5]) provides one equation involving these functions. So, we need 
another 2A^ — 1 constraints to fully determine W(x,a). We will consider the above ansatz for various values of 
N. Throughout this work, we shall use the following notation: 



(A + nh) 2 - A 2 . 
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• Lower case Greek letters will denote a- and x- independent constants; 

• Upper case Latin letters with a in parentheses will denote a-dependent, but x-independent "constants." 



N—l - Single- Term Superpotentials: 

W(x, a) = A(a) X(x) 
Substituting this ansatz into Eq.© leads to 



, x dA(a)\ „ .. . dX(x) dq(a) 
2 (A(a) X(x)) X{x) - 2A(a) , v ' + = 



da J dx da 



j2 dA dX\ dg 



Suppressing the x and a dependence : 2 A { X — — H — — = (13) 

da dx J da 



H(a) 



In Eq. tfTB")) . we see that the expression (A" 2 ij£ — can at most be a function of the parameter a, and hence, 
we set it equal to H(a). For H(a) ^ this leads to 



da 



Since, X(x) is independent of a, we must have, H (a) ~ d ^ a - ) and H (a) ~ ( d . These two conditions can 
only be met if H (a) and d "^ a ^ are constants, independent of the parameter a. Thus A must be a linear function 
of a: A — /ia + /3, and H (a) = 7 where the Greek-letter constants are all et-independent. 

This leads to the superpotential W(x, a) = — (/ia + /3)y^"tanh (^7 /i (x — xo)) . By a scaling (7/1 = 1) and 
a shift (xq = 0) of x, we can write giving 

W(x,a) = A tanha; , (15) 

where A = ^ — aj . This superpotential represents unbroken SUSY and holds finite number of bound state. It 
is a special case of the Scarf II and the Rosen-Morse superpotentials p] [2] • Note that, for negative values of H, 
instead of the hyperbolic trigonometric function, we get the superpotential W(x, a) — a tan x . Note also that this 
formulation includes terms like cot and coth as solutions of the same differential equation for different values of 
H(a) and ^p. For A — constant, Eq. (fT3|) generates 

W(x,a) = ^x , (16) 

the superpotential for the harmonic oscillator, where we have set H (a) = — ^lo. 

This exhausts the solutions for the N = 1 ansatz. Note that by definition it can only give one of the terms in 
any superpotential. 



6 For H (a) = 0, the solution of the above differential equation is X ~ l/x which is in fact a case of broken SUSY, and holds no bound states. 
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N—2 - Two- Term Superpotentials: 



W(x,a) = Y^Ma)Xi{x) 

»=i 

= Ai (a) X x (x) + A 2 (a) X 2 (x) 

We would like to very strongly emphasize that neither the functions Ai(a) and A 2 (a) , nor X\(x) and X 2 {x), 
should be linearly proportional to each other, otherwise this case will reduce to the case of N — 1 considered above. 

Substituting the above ansatz in Eq.© leads to 

2 (A Xl + A, X,) (*, + « ^ + = U 

Expanding it, we get: 



da da / da \ dx dx 



d -A-i v2 dXi dA 2 2 n a dx z , o v v { a d A 2 , A dAi\ t dg(a) 



2Ax-r± X^-2Ax-r ± +2M - r =-X^-2Aa-r ± +2XiX 2 A-t^+^-t^ + =0 . (17) 

v da ^ ax ^ da ^ dx y \ da dc 

Term#l Term#2 Term#3 Term#4 Tcrm#5 

This is of the form 

6 

i 

where 

F 1 =2A 1 ^; G 1 = (X 1 ) 2 

F 2 = -2Av, G 2 ^^ 



F 3 =2A 2 ^; G 3 = (X 2 ) 2 

F 4 = -2A 2 ; G^^ 

^=2(A^+^2^); G 5 ^X,X 2 

F 6 =dg/da; G 6 = l 



(18) 



This linear combination of various functions of x must yield us functions X\ and X 2 that are independent 
of a. This severely constrains the *4i's and Xi's. In fact, in the appendix we prove the following theorem: If a 
combination of FiGi is irreducibly constant, i.e., no smaller combination is independent of x, then all Fi's must 
be proportional to each other. We will be making extensive use of this property throughout this paper, where we 
utilize the above constraints to obtain all two-term superpotentials. 

In the notation we will be using, we will refer to Ft Gi as the i-th term; e.g, 2Ai X± as term jf 1, 2A\^-^- 
as term # 2, etc., as explicitly shown in Eq. (fTTj) . Of these terms, the sum of the first five terms is ^-independent 
since the sixth term, clearly does not depend on x. However, a subcollection of the first five terms may also 

add up to a constant. If, for example, terms 2, 3 and 4 add up to a constant and no smaller group of them is a 
constant, we would call such a group irreducibly constant and denote it by {2,3,4}. 
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Our general method will be driven by the following constraints: 

• Any irreducible term or set of terms that is irreducible must be at most only a-dependent. 

• In any irreducibly constant subgroup, the functions of Ai must be proportional to each other, as proved in 
the appendix. 

• For those cases where only a single term is irreducible, then either the function of Xi must be a constant, or 
if it is not, then the function of Ai must be zero, to preserve the dependence on a alone. 

1. Let us first consider the case that the sum Fi(a)Gi(x) is irreducible; that is, no smaller group of terms 
is constant. We will denote this group by {1,2,3,4,5}. In this case all F% are proportional to each other, as 
proved in the appendix. The proportionality of F 2 and F4 implies that the Ai's are proportional and hence 
this reduces to the N = 1 case. We would like to emphasize this feature. Whenever an irreducible group 
contains both term #2 and term #4, we have Ai ~ A 2 , and hence a case of N = 1. 

2. Now let us consider the case where there are four-term combinations independent of x. This also implies 
that there is one ir- independent term, which leads to severe constraint on X^s. The remaining 4-term group 
has all corresponding F^s proportional to each other. 

There are five such possibilities, of which, due to symmetry between the first and third terms and between 
the second and fourth, two of the terms are similar to others. For example, the condition U X\ =constant" is 
equivalent to "X 2 =constant." 
This leads to three cases: 

(a) {2,3,4,5} + {1} 

„ „ dXi „ „ dA> „n „ , dX 2 „ „ „ / „ dA 2 . dAi \ . 
- 2Ax—r- + 2-^2 —r 1 - X 2 ~ 2^2— + 2 X X X 2 Ai—^ + A 2 — L is an a-dependent constant. 
dx da dx \ da da J 

Since terms #2 and #4 appear in the same irreducible combination, the problem reduces to the N = 1 

case. 

(b) {1,3,4,5} + {2} 

d«4l ^2 r, A dA 2 ^ 2 n A , r, V V { A ^ A 2 ( rf^l 



2Ai -r^ X-f + 2A 2 —r^ X 2 Z - 2A 2 -r^ + 2 X X X 2 Ai-r=- + A 2 - r L = fli(o) (19) 
da da dx \ da da ' 

and [2-4i^] = - - H x (a). From the irredicibility of {1,3,4,5}, we have: 

. dAi . dA 2 , ( . dA 2 , dAi\ TT , , 
Ax —± ~ A 2 ~ A 2 ~ A 1 —^- + A 2 - r J - ~ H{a) . 



da da \ da da 

From A 2 „ A 2 , we get A 2 = aa + (3. From A x ^ ~ A 2 ; W e get (A) 2 = (A 2 f + S = 
(aa + l3) 2 +5. However, from +^2^) ~ A 2 = {aa+(3), we get A x = 7+Ccxp (-a 2 /2 - {pa)/ a) 

These two expressions for A\ are consistent only if it is a constant independent of a. But (A 2 ) 2 = 
(Ai) 2 — 8, so A 2 is constant as well. Constancy of both A\ and A 2 is a special case of A\ ~ A 2 . Thus, 
this problem reduces to the N = 1 case. 



7 As we did in the N = 1 case, we must distinguish between constants that depend on a, and those that are a- and ^-independent constants. 



7 



(c) {1,2,3,4} + {5} 

2Ai *± X? -2A^+ 2A 2 X* - 2A 2 ^ = H(a) , (20) 
da dx da dx 

and [2X 1 X 2 (Ai^£ L + A 2 ^)] =- - H(a). From the second and fourth terms of Eq. QHft, the 

irreducibility of {1,2,3,4} implies that A\ ~ A 2 . Thus problem reduces to the N = 1 case as well. 

3. There are ten possible ways to set three-term groups equal to a constant. However, due to the symmetry of 
terms with respect to X\ and X 2l we see that there are only six combinations we need to consider. We list 
them all below. They are: 

(a) Terms {1,2,3} add up to an a-dependent constant. 
Let us first consider {l,2,3}+{4,5}. 

2A X ^ X, 2 -2A^ + 2A 2 ^ X 2 2 = H(a). (21) 
da dx da 

Since we have Ai ~ A\ ~ A 2 , we find from the first condition that A\ must be a linear function 

of a, i.e., Ai — aa + (3, and from the second condition that (A 2 ) 2 = 7 {Ai) 2 + 6 = 7(0? a + [3) 2 + 6. From 

the constancy of {4,5}, we have 

_ 2A 2 ^ +2X 1 X 2 (A 1 ^+A 2 ^-)=J(a) . (22) 
dx \ da da J 

This gives A 2 ~ (Ai^ + A 2 ^) = k (Ai^ + A 2 ^) . Substituting A x = (aa+(3), we get the fol- 
lowing differential equation: k (aa + (3) + (k a — 1) A 2 = 0, which is solved by A 2 — (aa + /3)^ _1 . 
Thus, we have two expressions for A 2 , both of which must be valid for all values of a. Their compatibility 
implies — = 2, 7 = 1 and 5 = 0. Then A 2 = A\. This then is an N = 1 case. 

Now let us look at the case of {l,2,3}+{4}+{5}. 

From, {4}, we get X 2 = a 2 x + (3 2 . Since the fifth term is constant, we have two choices: it is either equal 
to zero or it is not. 

• Let us first consider the case that the fifth term is zero. As in the previous case, from {1,2,3} we 
have Ai^ ~ Ai ~ M cL ^\ i.e., A\ = aa + (3, and that (A 2 ) 2 = 7 (Ai) 2 + 5. But since the fifth 
term is zero, we also have Ai ~ A^ 1 ■ These two conditions can only be met if both Ai and A 2 are 
constants independent of a, and the problem reduces to the N = 1 case. 

• Now let us assume the fifth term is a non-zero a-dependent constant. Then X\ ~ \jX 2 . We have 
{1,2,3} and {4} separately equal to a-dependent constants. There are two possibilities for the fourth 
term: it is either zero or a non-zero constant. In the first instance, = 0, that makes X 2 a 
constant. This implies that {1,2,3} is reducible, and hence will be considered later. However, if 
the fourth term is not zero then X 2 = ax + (3. This implies that X\ = a ^ + fj ■ Then the irreducible 
three-term {1,2,3} implies: 

2M ^ (-Q 2 - JUx-^U + 2A 2 d -^(ax + f3? = M(a). 
da \ax + (3J (ax + (3) da 
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This algebraic equation has a nontrivial solution for all values of x only if 2A 2 = 0- However, 
that reduces {1,2,3} to {1,2}+ {3} and hence will be considered later. 

(b) Terms {1,2,4} add up to an a-dependent constant. 
I.e., 

2A^X* - 2Aj^ - 2A 2 ^ = H(a). 
da dx ax 

Then from the irreducibility of {1,2,4}, A\ ~ A 2 : an N = 1 case. 

(c) Terms {1,2,5} add up to an a-dependent constant. 

, dAi 9 . dXi ( , dA 2 . dAi\ , , , , 

da ax \ da da J 

The remaining two terms could add up to a constant irreducibly or separately. Thus, there are two 
possibilities for this case: {l,2,5}+{3,4} or {l,2,5}+{3}+{4}. Before we analyze these cases, we note 
that from Eq. (H3J), we get Ai ~ Ax ^ ~ (Ai ^ + A 2 ^). Thus, we have Ai = a x a + fi x . 
Now, we first consider the case of {l,2,5}+{3,4}. From {3,4}, we have, 

A 2 ^Xi-A 2 ^=C(a). (24) 
da dx 

From Eq. (El, we get A 2 ^ ~ A 2 , which gives A 2 = a 2 a + j3 2 . Substituting expressions for A\ and 
A 2 in A\ — (M^jjr + A 2 ^jjfc*) = a 2 Ai+ ai A 2 implies that A 2 ~ Ai, thus reducing this to an TV = 1 
case. 

Next we consider the case of {l,2,5}+{3}+{4}. From {1,2,5}, we still have A\ = aia + Pi and A\^- ~ 
(A 1 ^~Ja L +"^ 2 ^3a 1 ) ~ Ai. = 7i A\. From {3} there are two possibilities. Either A 2 is constant or X 2 is 
constant. If A 2 is constant, from {1,2,5} we find that A\ is also constant, and hence we have an N = 1 
case. Now we consider the case X 2 = 72 is constant. The irreducible combination now leads to the 
differential equation 

2 dX x H{a) 

a ^ -^x~ + ^ Xi= za;> 

whose solution is ^tanhx — • To determine A 2 , we start with the equation (Ai^ 1 + A 2 ^ L ) = 
J1A1, i.e., 

dA 2 

(aia + /3i)— \-aiA 2 =7i(aia + /3i), (25) 

da 



or 



Eq. (US]) is solved by 



^ + ^hf-4 2 -7 1 = 0. (26) 
da a,\a + p\ 



A2 = 7i (|aiQ 2 + /?ia) +Ci ^ 
aia + pi 



Solving for the superpotential W(x, a) = A\X\ + A 2 X 2 , we get 
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W(x, a) = (aia + fti) tanhx + 
= A tanh a 



t - 111 

"»1 2qi 



a,\a + /?i 

1 ' 



where A = a±a + fli and _B = £1 — 2^ 



2qi • 

(d) Terms {1,3,4} add up to an a-dependent constant. This is the same as item (a). 

(e) Terms {1,3,5} add up to an a-dependent constant. 

d Al ,, 2 l0 /l V2 . r, « V ( A d M , , d M 



2A X X{ + 2A 2 X 2 + 2 <*i*2 A x —^ + A 2 - r ± = H(a). (28) 
da da \ da da 



Eq. m implies that we have Ai^ ~ ~ ( A ^+^ 2 ^). „ ^ 2 ^i g ive s 

(A) 2 = ai (A 2 f + (3i. The proportionality condition Ai^± ~ (Ai^ + A 2 ^) generates {Aif = 
a 2 (AiA 2 ) + (3 2 . If pi or j3 2 vanishes, we get A\ ~ A 2 . If fii and /?2 are different from zero, these two 
algebraic equations for Ai and A 2 determine them to be constants. Thus, in both situations, we have 
an N = 1 case, 
(f) Terms {1,4,5} add up to an a-dependent constant. 

Let us first consider the case {1,4,5} + {2,3} . Thus, we have 

2Ai X? - 2A 2 ^ + 2 Xl X 2 (Ai^ + A 2 ^) = H x {a). (29) 
da dx \ da da J 

and 

2A 2 ^ X 2 2 -2Ai^= ff a (o). (30) 
da dx 

Eqs. ([29|) and ([SOD imply A x *fe ~ At ~ {Ai^+A 2 ^), and yl 2 ^ ~ ^ respectively. Thus, 
we have ^ ~ £ and *fc ~ From ~ + A^), we get A 2 = 5 X A X + 

Substituting above expressions for ^j 1 , 4?> an d -^2 into -4 2 ~ (Ai^ 1 + -A^j; ), we get the following 
algebraic equation for A\: 

2 



S 2 _ ^ , + £ 



cMi + -r- = , • 

Ai + % Ai 

whose solution A\ equals a constant. However, this then makes {1,4,5} irreducible, hence will be 
considered later. 

Now consider the case {1,4,5} + {2} + {3}. From {3}, we have A 2 a constant because X 2 constant 
makes {1,4,5} reducible. Using this in {1,4,5}, we get A\ = 71 a + Si and A\ = 72 A\A 2 + S 2 . Since 
A 2 is constant, from the second condition we get Ai constant. The problem thus reduces to the N = 1 
case. 

(g) Terms {2,3,4} add up to an a-dependent constant. This is the same as (b). 
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(h) Terms {2,3,5} add up to an a-dependent constant. This is the same as (f). 

(i) Terms {2,4,5} add up to an a-dependent constant. As stated earlier, the presence of term #2 and term 
#4 in {2,4,5} implies A x ~ A2 and the problem is immediately reduced to an N = 1 case. 

(j) Terms {3,4,5} add up to an a-dependent constant. This is the same as (c). 

4. As we set out to now consider all irreducible two term cases, let us recognize again that term #1 in Eq. (fl~7|) 
is equivalent to term #3 and term #2 is equivalent to term #4. There are several possibilities. 

(a) . {l,2}+{3,4}+{5} 

{l,2}+{3,4}+({5} = 0) 

The first and second terms add to an a-dependent constant: 2Ai X x 2 — 2A x ^r- = H(a). The 
proportionality of the coefficients on the left hand side imply = a. Then A\ = aa + rj. Similarly, 
the constancy of the combination of the third and fourth terms gives A2 = 012a + t)% 
But from {5}, A2 ~ l/A x ; i.e., a^a + r]2 ~ l/(a x a + r) X ). For this to be true for all values of a, 
oci = ct2 = 0. Then both Ai and A2 are a-independent constants: a special case of A2 ~ Ai, and 
thus an N — 1 case. 

{l,2}+{3,4}+({5} ^ 0) In this case, instead of A 2 ~ 1/A X , from {5} we have X 2 ~ l/X x . The 
differential equations from {1,2} and {3,4} give respectively X x = tanhx and X% = cothrr. As a 
result, the superpotential for this case will be of the form: 

W(x, a) = (a x a + rj x ) tanh (7 x) + {ot2a + 772) coth (7 x) . (31) 

With the identification of (a x a + rji) = A and (a^a + — —B, we get Atanhr — B cothr. This is 
the Poschl-Teller II superpotential. This also produces, with different constants, the Poschl-Teller 
I superpotential [2], 
. {l,2}+{3,5}+{4} 

The constancy of the fourth term implies that ^j 2 - — a; i.e, X% = ax-\- (3. We cannot consider the 
case of a = here because that would make X2 constant, and thus {3,5} would be reducible to {3} 
and {5}. We will consider this later. From the irreducible {3,5}, we thus obtain 

A 2 ^ (X 2 ) 2 + X X X 2 4- {A X A 2 ) = D(a). (32) 
da da 

From this we get A 2 ^ ~ ^ {AiA 2 ), i.e., Ai = 0iA 2 + This leads to 

{X 2 f + 6 X X X X 2 = k. (33) 

For k — 0, we have an algebraic equation which reduces to X\ ~ X2; thus an A = 1 case 

For k ^ 0, from {1,2}, we have A x ^ {X x f ~ A x ^ = J (a), i.e., A x ^ ~ A x . This gives 

Ai = 7a + v, which yields 

/ v \2 dX x 
l{Xi) - — = A. 
dx 
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For A = 0, we get X\ = — x 2 Xa ■ This is not compatible with the condition given in Eq. (|33[) and 
X2 = ax + [3. Thus there is no solution. 

For A 7^ 0, we get two cases again. 7 could be equal to zero or a non-zero constant. For 7 = 0, 
{1,2} becomes reducible, and hence will be pursued later. For 7 7^ 0, the solution for X\ will be a 
function of the type tan, cot, tanh, or coth, none of which is compatible with X2 = olx + [3 and Eq. 
(|53"1) . Thus there is no solution. 

. {l,2}+{4,5}+{3} 

Constancy of the term {3} implies that either X2 = £, a constant or the derivative of A2 'zero, 
i.e., A 2 = a 2 . 

— If X2 is a constant, then {4,5} becomes reducible and will be discussed later. 

— If A2 is constant, from {4,5} we have A\ — aa + (3. From {1,2}, we have 

aX 2_dX 1 = x (34) 
dx 

For A = 0, we get X\ = — ax 1 _ ( ^ 1 Substituting in the equation (from {4,5}) 

dX 2 



dx 

and solving for X 2 , we get 

X 2 = 



aX l X 2 = v , (35) 



^(fz 2 -Ci*) +C2 



a x — Ci 

The superpotential is then given by 

aa + P+^ + (2 & . 

y 5-(ax-Ci) , 

a x — Qi la 

which for aa + (3 + +£2 = ^ + 1, — 2 = % an< ^ a x ~ Ci = r i becomes 

W{r,a) = -ujr , 

2 r 

the three-dimensional harmonic oscillator. 

For A 7^ 0, we have X\ is of the form tana;. Then solving for X 2 in Eq. (|35[). we get A2 of the 
form sec a;. Hence, the superpotential generated is of the form 

W(x, A) = A tan x + B sec x. 

Depending on values of constants a, v and A, Eqs. (l34|) and |35|) also produce a shape invariant 
superpotential: 

W(x, a) — A coth a: — -Bcosech x. 

This is the generalized Poschl-Teller superpotential pQ. Similarly, Eqs. (|34|) and (|35p can also be 
used to generate the Scarf I (A tan x — Bsec x) and Scarf II (A tanh x -I- i3sech x) superpotentials. 
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{l,2}+{3}+{4}+{5} 

From {1,2}, we have Ai ^ Xf - Ai ^ = M(a). This gives ^ = a, constant; thus, Ai = aa+/3. 
For M(a) ^ 0, we then get X\ <~ tanh or X\ ~ coth or their trigonometric counterparts; i.e., tan 
and cot. For M(a) — 0, the solution is X\ = — §■ 

For X2, let us turn to the remaining terms. The fifth term is either zero or a non-zero constant. We 
will consider them separately. 
{l,2}+{3}+{4}+({5} = 0) 

From {3}, either X 2 or A2 is a constant. If X 2 is constant, with A\ <~ l/A 2 from {5} =0 and 
M{a) 7^ 0, we generate the following superpotentials: Rosen Morse I (— -j — ^4cota;), Rosen Morse 
II (f + Atanha;), Eckart (f - A coth x) . For M(o) = 0, we get instead W(x,a) = - a ( aa +^ + 
a a + p - With the identifications x — > r, a(aa + /?)—> (Z + 1) and = 2 (f+ 1) > we obtain the 

Coulomb superpotential. 
{l,2}+{3}+{4}+({5} + 0) 

From ({5} 7^ 0), we have X\ <~ l/A^. Thus, A2 = constant is no longer a solution. From {3} and {4}, 
we find that X 2 must be a linear function of x and A 2 must be a constant. Thus, only solution for X\ 
that is compatible with X\ ~ 1/ A" 2 and linearity of A2 is <%i = — ^ . Thus the superpotential generated 
is W(x, a) = - a( °° +(i) + 72;. With the identification of x -> r, a(aa + /?) -> (Z + 1), and 7 -> ± w, 
we get the superpotential of the three-dimensional harmonic oscillator: 

W{r,a) = - ujr — , 

the same as {l,2}+{4,5}+{3}. 
{l,3}+{2,4}+{5} 

From {2,4}, this combination immediately reduces to an N = 1 case. 
{l,3}+{2}+{4}+{5}. 

Let us first consider {l,3}+{2}+{4}+({5} = 0). 

Since {2} and {4} are constants, we have X\ = ot\ x + (3\ and X 2 = a 2 x + f3 2 . From {1,3}, we have 

(ai x + Pxf + n{a 2 x + (3 2 f = S. 

The above equation is valid for all values of x, only if a\ x + (3\ = \J—k (a 2 x + f3 2 ); i.e., X\ <~ X 2 . 
Thus, we have a case of N = 1. 

This satisfies {1,3} only if X\ and X 2 are proportional. However, that again reduces the problem to 
the N = 1 case. 

{l,3}+{2}+{4}+({5} ^ 0), 

In this case since and X\ X 2 are all constant, we get we get the trivial case of constant 

functions X\ and X 2 . 
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{l,3}+{2,5}+{4} 

From {4}, we have X 2 — -fx + 5. Since terms 1 and 3 add irreducibly to a constant, we have 
(Ai) 2 + k(X 2 ) 2 = v. From this, we get X\ = \jv — k (72: + 5) 2 . Substituting these expressions of 
X\ and X 2 in the combination {2,5}, we see that it is satisfied for a range of value of x only if 7 = 0. 
Thus, we get the trivial solution of X 2 and X 2 both constant. 

{l,3}+{4,5}+{2} is equivalent to the case {l,3}+{2,5}+{4} considered above; viz, the trivial solu- 
tion. 

{l,3}+{2}+{4}+{5} 

First consider {l,3}+{2}+{4}+({5}=0). Since the sum of the first and third terms, 2Ai X 2 ^ + 
2A 2 X 2 ^ = B(a), we have A\ = eA\ + 7. But from ({5}=0) we have A 2 ~ 1/Ai. Thus the only 
solution is that both Ai and A 2 are a-independent constants, hence the special case of A 2 ~ A\; 
viz, the N = 1 case. 

Next we consider {l,3}+{2}+{4}+({5}^0). 

From ({5}^0), we have X 2 ~ 1/X\. Terms {2} and {4} imply that X\ and X 2 are linear functions of 
x. Compatibility among {2}, {4}, and {5} thus requires that X\ and X 2 both be constants, therefore 
we get the trivial solution again. 
{l,4}+{2,3}+{5} 

Let us first consider {l,4}+{2,3}+({5} = 0). In this case, from {1,4} we have ~ and from 

({5} = 0) we have Ai <~ l/A 2 . The product of these two conditions give A\ ~ constant, i.e., 

d< "~^a — Q: ' a constant. This gives Ai = (aa + fi) 1 ^ 3 and A 2 = (aa + /3)~ 1//3 . However, from {2,3}, 
we also have ~ A\. This last constraint implies that we must have a = 0, i.e., A\ and A 2 are 

constants; viz, the N = 1 case. 
Now let us consider {l,4}+{2,3}+({5} ^ 0). 

From ({5} ^ 0), we have X x = j/X 2 . From {2,3} we get ^ + aiAf = (i u and from {1,4} we have 
^ + a 2 Xl = (3 2 . Differentiating X 1 = j/X 2 , we get ^ = -7^/ (X 2 f . Eliminating X x and its 
derivative from these two equations, we get an algebraic equation for X 2 , whose solution is X 2 = 5. 
This, along with X\ = j/X 2 , reduces the problem to a trivial case of X\ and X 2 equal to constants. 
{l,4}+{2,5}+{3} 

Constancy of the term {3} implies that either X 2 is constant or = 0. 

— First consider the case that X 2 is constant. This gives {4} = 0; i.e., {1,4} becomes a reducible 
case to be considered later. 

— For the other possibility, = 0, A 2 is a constant. Then from {1,4}, wc get 4- (A\) 2 = Q\. 
Which gives A\ = \/0\a + 9 2 . However, from {2,5} we get Ai ~ j^Al, i- e -7 -^l = 7ie 72Q + 73- 
These two expressions are not compatible unless both Ai and A 2 are constants, in which case 
N = 1. 

{l,4}+{3,5}+{2} 
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{2} implies that ^± = constant. So, X\ = a\x + /3\. 

- ai = 0. 

Thus X\ =/3i. This makes {1,4} reducible. 

- ai + 0. 

From {1,4}, we have X 2 ~ (aix + /3i) 3 + 7. 

But {3,5} implies A^ 2 <~ A4A2 + S. This expression is not compatible with the values of Xi and X 2 
above, unless ai — 0; i.e., Xi and X2 are constant, the trivial case once again. 
. {l,4}+{2}+{3}+{5} 

As we saw above, constancy of the term {3} implies that either Xi is constant or Ai is constant, and 
Xi constant will not be considered, since it makes {1,4} reducible. Thus Ai is constant. For {5} = 0, 
Ai ~ 1/Ai = a constant, thus reducing to an N = 1 case. For {5} 7^ 0, we have X\ <~ l/Xi. Term 
{2} constant means that X x = ctix + f3i, hence X 2 <~ l/(a 1 x + (3^. But from {1,4}, = eX i + 
So Xi — e(aia; + /3i) 3 /(3Q;i) + (5x. These are incompatible expressions; therefore, there is no solution. 

• {l,5}+{2,3}+{4}. This is equivalent to {1,4} + {3,5} + {2} analyzed earlier. 

• {l,5}+{2,4}+{3}. The presence of {2,4} immediately implies an N = 1 case. 

• {l,5}+{3,4}+{2} is equivalent to {1,2} + {3,5} + {4} analyzed earlier. 

• {l,5}+{2}+{3}+{4} From {2} and {4}, we have X x = a Y x + fii and X 2 = a 2 x + (3 2 . Then from 
{3}, we have two possibilities: either X 2 is a constant or Ai is constant. 

— We first consider the case Ai constant. This implies, from {1,5}, Ai is constant as well. This is 
an iV = 1 case. 

— Now we consider the case Xi = (3i, a constant. Since we have X\ = ol\x + fli, from {1,5} we get 

dAl . n \1 n I n \ AA\Al T/ N 

1 ~cto~ ^ + ^' + ^ 2 ^ + ^ da = ^' 
This equation can be satisfied for all values of x, only if ai = , in which case X\ is also constant; 
viz., a trivial solution, or if Ai and Ai are independent of a, reducing this to an N — 1 case. 

• {2,4}+{l}+{3}+{5}. {2,4} gives Ai - A 2 , an TV = 1 case. 
. {2,5}+{l}+{3}+{4} 

Since term {1} is constant, we have two choices. It is either equal to zero or a non-zero constant. 

• If it is non-zero, that would imply that X\ is a constant. This reduces {2,5} to {2} +{5} and will 
be considered later. 

• On the other hand, if term {1} is zero, we must have Ai a constant (since X\ — is not acceptable.) 
So we choose Ai = r/. From {2,5} we then get = a, so Ai = aa + j3. Further, from {3} we 
have X 2 = 5, a constant , since Ai is not. This then gives (4|r ~ 018X1) = fi. The solution of this 
equation is given by X\ — Ke^ aSx ^> — ^ . Thus the resulting superpotential W(x, a) = AiX\ + -42 A2 
is given by 

W(x,A) = A-Be~ ux , (36) 
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where A = {oca + (3) — B = rj, v — aS. For a < 0, this gives the Morse superpotential. 

(g) {l}+{2}+{3}+{4}+{5}. 

Let us consider various possibilities. 

({1} = 0) implies Ai is a constant. ({1}^0) implies X\ is a constant. ({3} = 0) implies A2 is a constant. 
({3}t^0) implies X2 is a constant. So, if both ({1} = 0) and ({3} = 0), we have A\ ~ A\, hence N = 1. 
Similarly, if both ({1}^0) and ({3}^0), we have X\ and X 2 equal constants, a trivial case. Thus, the 
only case that we need to consider is ({1}^0) and ({3} = 0). This equivalent to ({1}=0) and ({3} ^ 0). 
In this case, we have A\ = v\ and X2 = $2- Then {5} yields v\f3iX\^. If ({5} = 0), we have A2 = 1^2] 
i.e., Ai ~ A2, an N = 1 case. If ({5} 7^ 0), X\ — f3\ and hence a trivial case where both X\ and X2 are 
constants. 

This exhausts all the possibilities for the N = 2 ansatz. 

In particular, we have now obtained all known translational shape invariant superpotentials. We have 
found no new ones, using the N = 1 and N — 2 ansatz. Table I lists these superpotentials, as well as 
one of the combinations of irreducible terms from which they were obtained @. 



Name 


superpotential 


N 


Combination 


Harmonic Oscillator 




1 




Coulomb 


e 2 l+l 
2(J+1) r 


2 


{l,2}+{3}+{4}+({5} = 0) 


Three dim. harmonic oscillator 


2 r 


2 


{l,2}+{3}+{4}+({5}^0) 


Morse 


A - Ber x 


2 


{2,5}+{l}+{3}+{4} 


Poschl-Teller I 


A tan r — B cot r 


2 


{l,2}+{3,4}+({5} * 0) 


Poschl-Teller II 


A tanh r — B coth r 


2 


{l,2}+{3,4}+({5} 0) 


Rosen-Morse I 


-Acotx- -f 


2 


{l ) 2}+{3}+{4}+({5} = 0) 


Rosen-Morse II 


Atanhx + -f 


2 


{l,2}+{3}+{4}+({5} = 0) 


Eckart 


-A coth x + -§ 


2 


{l,2}+{3}+{4}+({5} = 0) 


Scarf I 


A tan x + Bsec x 


2 


{l,2}+{4,5}+{3} 


Scarf II 


A tanh x + _Bsech x 


2 


{l,2}+{4,5}+{3} 


Generalized Poschl-Teller 


A coth x — _Bcosech x 


2 


{l,2}+{4,5}+{3} 



Table 1 

As we have seen, there are multiple combinations of irreducible terms which obtain the same superpotential. 
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Conclusions 



Until now, finding shape invariant superpotentials was left to excellent intuitions or trial and error. The 
procedure described here is an ab initio method for generating superpotentials that are shape invariant, and 
hence exactly solvable. 

We have transformed the translational shape invariance condition into a nonlinear partial differential equation. 
With this step, we have transformed a non-local difference-differential equation into a differential equation which 
must be satisfied by all translational shape invariant superpotentials. We have also constructed a variant of the 
"separation of variables" method to find solutions of this nonlinear partial differential equation, and have found 
that the list of solutions generated includes all shape invariant superpotentials, as given in Ref.[5]- 

This work can be extended in several directions. On the one hand, other forms of shape invariance can be 
investigated, to ascertain whether a similar transformation of difference equation into differential equation is 
possible. In particular, while we have not considered the case of multiplicative shape invariance, we believe that 
this work can be extended in that direction. 

On the other hand, in an earlier work [5] we argued, based on group theoretical considerations, that all of the 
known cases of translational shape invariance are the only ones which can exist. We have verified this claim for 
N = 1,2. However, it remains to be checked whether our claim holds for larger values of N. 
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Appendix 

Linear dependence of coefficients Fi for an irreducible set 

The nonlinear partial differential equation that we want to solve is given by 

2 Al X? ~ + ^ ^ X} - 2^ +2^ 2 ( A £± + A 2 ^-) = - ^ (37) 




Term#l Tcrm#2 Tcrm#3 Term#4 Term#5 Tcrm#6 

This is of the form 



J2Fi{a)Gi{x) = (38) 



where 



F 1 =2A 1 ^; G l = {X l f 

F 2 =-2A i; G^ d -£ 
dA 2 . n — I v \ 2 



F 3 =2„4 2 ^; G 3 = (A- 2 ) 

F^~2A 2 ; 

F 6 =dg/da; G 6 = 1 



(39) 



Eq. (f3"8"|) shows that not all Fi are linearly independent. In particular, it states that at most five Fi can be 
linearly independent and the dimensionality of the space they would span could at most be five. We will show 
that if the expression irreducibly adds up to a constant, the space spanned by Fi would have the dimensionality 
of one; i.e., all F$ must be proportional to each other. 

Let us first consider an expression consisting of just two elements: {FiGi, -F2G2}, that add up to a constant 
irreducibly; i.e., Fi(a)G\(x) + F2(a)G 2 (x) = H(a). The irreducibility for this case implies that Fi(a)Gi(x) and 
F2(a)G2(x) cannot be separately constant. This implies that neither G\(x) nor Ga(a;) can be ^-independent 
constants, and F\(a) and -F 2 (a) cannot be equal to zero, otherwise the two terms would be reducible to one. 

Dividing the above equality by H(a) and defining T\{a) = §^ ■ •^" 2 ( a ) = (a) 1 

Ti(a)Gi(x) +J r 2(a)G 2 {x) = 1. (40) 

Since x and a are real variables, the above expression must be valid for infinitely many values of both x and a. 
Let us consider a\ and 02, two arbitrarily chosen values of a. For them, we have 

^i(ai)Gi(a;)+^ 2 (oi)G 2 (a:) = 1 

Fl^G^x) + F 2 {a2)G 2 (x) = 1. (41) 

Imagine a plane where G\ {x) and G2 (x) are x and y coordinates respectively. Each of the equalities expressed 
in Eq. (|4"Tj) are represented by a line on this G1-G2 plane. If these two lines intersect, we will have a solution 
for both G\{x) and G 2 (a;); in other words, each would be a determined constant. However, our hypothesis of 
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irreducibility states that we should not be able to determine such solution; i.e. Eq. (|41[) should be invertible 
This implies that we must have the lines parallel; i.e., both must have the same slope: 

•^(ai) ^2(02) 



J 7 ! (ax) Fi{a 2 ) 



(42) 



In other words, the ratio is independent of the argument a, and hence must be equal to an a-independent 

constant. This implies F\ and F 2 are proportional to H. 

Let us now consider an example that consists of three terms that irreducibly add up to a constant. So, wc 
have 

^ 1 (a)G 1 (x) + ^ 2 {a)G 2 (x)+^ 3 (a)G 3 (x) = 1. (43) 

Again, we can define J-i(a) = for i = 1, 2, 3, and obtain the following set of three equations for three arbitrary 
values of a: 

^i(a 1 )G 1 (x)+^2(ai)G 2 (x)+^ 3 {ai)G 3 (x) = 1 
Fi{a 2 )G x {x) + F 2 {a 2 )G 2 {x) + T 3 {a 2 )G 3 {x) = 1 

TMG^x) + T 2 (a 3 )G 2 (x) + T 3 (a 3 )G 3 (x) = 1. (44) 

In this case, consider a three dimensional space where Gi(x) (i — 1,2,3) are along coordinate axes. Eqs. (|4~4"|) 
represent three planes. If these planes intersected, the functions G{(x) will get determined and our expres- 
sion would reduce to two-term or single-term expression. Since we stipulated that our three-term expression 
was irreducible, these planes must all be parallel. This implies that the direction ratios of their normals; i.e., 
{Ti(ai), T 2 {ai)i T 3 (a*)} must be proportional: 

•FjXgi) _ ^(ai) _ •7g(ffli) 
Ti(a 2 ) T 2 (a 2 ) T 3 (a 2 ) 

must all be equal for different values of a,; . This equation can be rearranged to yield 



(45) 



This means that ratios such as ^fej do not depend on a. Thus Ti{a) ~ J- 2 (a) ~ ^3 (a). This implies F l5 _F 2 , 
and i 7 ^ are proportional to i7. 

The method we used for two terms is identical to the one we used for the three-term case and can therefore be 
similarly extended to higher dimensional spaces for any of the four- or five-term irreducible expressions we have 
used in this paper. 



19 



